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Abstract 

The null-brane space-time provides a simple model of a big crunch/big bang sin- 
gularity. A non-perturbative definition of M-theory on this space-time was recently 
provided using matrix theory. We derive the fermion couplings for this matrix model 
and study the leading quantum effects. These effects include particle production and 
a time-dependent potential. Our results suggest that as the null-brane develops a big 
crunch singularity, the usual notion of space-time is replaced by an interacting gluon 
phase. This gluon phase appears to constitute the end of our conventional picture of 
space and time. 



^ email address: 
^ email address: 
^email address: 



ej m@theory. uchicago . edu 
robbins@ theory, uchicago .edu 
sethi@theory.uchicago.edu 



1 Introduction 



One of the major gaps in our understanding of string theory is how it resolves space-hke 
singularities. A related issue is how to formulate dynamics in time-dependent backgrounds 
non-perturbatively. The various gauge/gravity correspondences subsume a resolution of 
space-like singularities when they are localized inside black hole event horizons. Implicitly, 
there must be encoded the experience of a freely falling observer who sees a time-dependent 
background culminating in a space-like singularity. One might then wonder whether there 
is a regularization of cosmological singularities via gauge theory. 

Recently, some examples of such regularizations via gauge theory have been concretely 
described in [1-14]. These holographic models, based on matrix theory [15], appear to 
contain the right degrees of freedom to describe physics at large blue shifts. In each case 
the semi-classical picture of the space-time singularity is replaced by a gluon phase. We 
will find further evidence for this picture in this work. Quantum effects in such models are 
quite fascinating have been examined recently in [16,17]. 

Our goal is to develop this matrix approach when applied to M-theory in the particular 
background known as the nuU-brane space-time [18,19]. The matrix model for this space- 
time was introduced recently in [20]. The nuU-brane is constructed using a particularly 
simple way of generating time dependence: namely, by orbifolding Minkowski space by a 
discrete subgroup of the Poincare group. This can yield a time-dependent background. For 
appropriate choices of orbifold group generators, this background can be singular. 

The nuU-brane is actually a non-singular solution of M-theory defined by choosing an 
R^'^ subspace of R^'-'^^ with coordinates 
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and the usual flat metric ds^ = —2dx'^dx + dx^ + dz^. We act on these coordinates by an 
element of the 4-dimcnsional Poincare group 
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When L goes to zero, the geometry has a null singularity at = 0. This space-time is 
depicted in figure HJ 




Matrix quantum meclianics 

Figure 1: The circle radius shrinks to a minimum size L at = 0. 



For L large compared to the Planck scale, this is simply a smooth time-dependent 
background of M-theory or string theory with a well-defined S-matrix. String perturbation 
theory on this background has been studied in [21-23] while the semi-classical stability of 
this space-time has been examined in [24, 25] . Not surprisingly, string perturbation theory 
breaks down as L ^ 0. Holographic descriptions of string cosmologies, in the spirit of 
AdS/CFT, have been studied in [26-28]. In the case of the nuU-brane, these descriptions 
involve space-time non-commutative Yang-Mills theories. 

By contrast, the matrix description is obtained by performing an additional light-like 
compactification x~ ~ x~ + 2ttR which commutes with the orbifold identification (0). This 
results in a discrete light-cone quantized description of the theory. In close analogy to the 
fiat space case [29], there is a decoupling limit [20] that reduces the matrix model to the 
low-energy dynamics of DO-branes moving in the orbifold Q. 

The matrix theory description can also be thought of as viewing the system from the 
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perspective of a highly boosted frame. It does not directly describe the space-time vacuum, 
but rather describes objects in space-time that carry large p"*", and therefore preserve at 
most one-half of the supersymmetries. In the case of the null-brane, the one-half of the 
supersymmetries preserved by the orbifold identification are incompatible with the one- 
half of the supersymmetries preserved by boosting along the light-cone direction. Thus 
the null-brane matrix model preserves no supersymmetry. We will see this explicitly when 
we evaluate the one-loop effective potential of the model, and compute the amplitude for 
particle production. 

We can understand the dynamics of a DO-brane in this orbifold background from the 
perspective of the covering space of the orbifold group action Q. On the covering space, we 
see a DO-brane together with an array of boosted images constrained to move in a manner 
invariant under the identification ((21). The interactions between the image branes arise from 
the usual velocity-dependent forces between DO-branes. 

Semi-classically, when the brane and enough of its images are contained within their 
Schwarzschild radius, a black hole forms. For finite L, the black hole eventually decays and 
we can compute a conventional S-matrix. However, as L becomes small, this black hole 
becomes larger eventually filling all of space [25]. We will find a quite different picture of 
this process from the matrix model which, in particular, includes degrees of freedom not 
seen in the semi-classical picture. 

T-dualization of this system of DO-branes along the z direction circle yields a Dl-brane 
matrix string theory description. The evolution from matrix quantum mechanics to matrix 
strings near the singularity is also depicted in figure H Physics near the singularity is 
controlled by 1 + 1-dimensional field theory dynamics rather than quantum mechanics. In 
principle, the future phase is a return to matrix quantum mechanics as the z circle becomes 
large again (and the T-dual field theory circle becomes small). 

However, strong quantum effects near the singularity can drastically alter the future 
description. Determining these quantum effects at one-loop in the matrix model is our 
primary task. We will find that a potential is generated in the matrix model which turns 
off at early and late times. This supports the consistency of the model since we need exact 
fiat directions in the far past and future to recover a semi-classical picture of space-time. 

The potential tends to attract gravitons more and more strongly as they approach 
x"*" = and as we tune L — 0. This effective potential as a function of the impact parameter 
for two gravitons is displayed in figure [21 As the gravitons approach each other, the off- 
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diagonal matrix degrees of freedom become lighter and our notion of space-time becomes 
fuzzy and is eventually replaced by matrix dynamics. At this level of approximation, it 
appears that there is no escape from this gluon phase as L — 0. 

Potential 

' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' — Impact parameter 
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Figure 2: The potential given in (|<S8|) as a function of impact parameter with various choices of L. 
The smaUer L, the deeper the potential well. The choice of parameters is described in section im 

The attraction of gravitons as and L go to zero is a reflection of their interaction with 
one another's images under the orbifold group. In fact, even a single graviton interacts with 
its own images and can generate a spacelike singularity at the nuU-brane neck at = 0; 
even if the vacuum does not evolve to a singularity, excited states can do so. The nuU-brane 
matrix model describes objects (such as gravitons) carrying longitudinal momentum in the 
null-brane background, rather than being a description of the nuU-brane vacuum directly. 
Thus even the ground state of the matrix model can evolve to a singularity when L is small 
enough. This property is related to the lack of supersymmetry of the null-brane matrix 
model. 

The formation of such singularities by excitations on the null-brane was analyzed by 
Horowitz and Polchinski from the point of view of the dual gravitational theory in [25]. 
They argued that a spacelike singularity will form in the null-brane orbifold whenever the 
gravitational radius of excitations exceeds the proper size L of the null-brane neck at a;"*" = 0. 
For d > A noncompact spatial dimensions, the singularity that results is a finite-sized black 
hole;^ outside the black hole there is no big crunch singularity. The gravitational radius of 
^More precisely, a black string which fills the z direction. After the neck is passed and the size of the z 
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the black hole is 



L2 ) 



(3) 



in terms of the effective d + 1 dimensional Newton constant Gn and momentum p"*" of the 
initial excitation. In particular, the size of the black hole grows to infinity as L — > 0. For 
d — 4:, no black hole forms for large L, and an infinite mass black hole forms for small L. 
For d < 4, a spacelike singularity of infinite extent always forms. 

The interesting regime of d < 4 is outside the reach of matrix theory. Matrix theory 
constructions describe gravity in terms of a dual field theory for d > 5, and a dual "little 
string theory" ior d — 5; there are no constructions for d < 5 (because the growth of the 
density of states is too rapid to be described by a field theory or nongravitational string 
theory). Thus, for finite L, matrix theory is only capable of describing the regime where 
the outcome of objects faUing through the neck of the nuU-brane is a finite-sized black hole. 
The size of that black hole grows to infinity as L/£p goes to zero. One might view the 
matrix theory as a description of the cosmological horizon that forms in this limit. 

Black holes are well-approximated in matrix theory at finite N provided their entropy is 
less than [31-33]. Therefore, the growing size of the final state black hole with decreasing 
L puts a lower bound on N 



in order to have the matrix model approximate well the bulk gravitational physics. A 

conservative order of limits sends N ^ oo first, then L — > to achieve a cosmological 

singularity. The dynamics may sensitively depend on this order of limits, however an 

intriguing possibility is that one could obtain a useful picture of the cosmological singularity 

by taking L ^ for finite and then subsequently the large N limit. The L ^ limit 

of the finite theory accesses the perturbative regime of the gauge theory; if this order 

of limits is germane, the cosmological singularity devolves to a kind of free gluon phase. 

If on the other hand the order of limits docs not commute, one will need to contend with 

the strong coupling, large physics of the gauge theory; still, the matrix model seems to 

provide a resolution of the cosmological singularity. 

This picture of the null brane is very much in accord with the model studied in [1, 17] 

in which time begins in a free gluon phase. While our understanding of the gluon phase 

circle starts to grow, a Gregory-Laflamme instability sets in [30] and this black string breaks up into one 
or more black holes. 




(4) 
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(free or otherwise) in the nuU-brane is still quite preliminary with much to be understood, 
it is exciting that the matrix model appears to be a complete description containing all the 
ingredients necessary for understanding the fate of the S-matrix and, indeed, the fate of 
space and time. 

2 The Matrix Action 

Our first task is to determine the complete action for the matrix description of M-theory 
on the nuU-brane. The bosonic couplings were determined in [20], but we also require the 
fermion couplings to compute quantum corrections. 

2.1 Parameters, scales and sectors 

It is important to keep track of the scales appearing in this system. The space-time physics 
depends on the matrix theory parameters {£p, R) as well as the orbifold parameters (L, A). 
The former are related to string theory parameters {is,gs) via 

ip = gl'Hs, R = gds. (5) 

We will express the matrix action in terms of string parameters and set = 1 for simplicity. 
The bosonic couplings in the matrix model describing M-theory on the nuU-brane with 
p"*" = N/R are given by a 1 + 1-dimensional U{N) gauge theory with action [20] 

S = J j drdaTi (^{Dox'y + {Doxf - L^^D^x'f - L^D^xf - V2\{x - tDox)F 

+ {L' + lx'r')F' + l[x\ x^] ' + {[x, x^] + ^^i^^O^) , (6) 

where < cr < 27t. Here r can be identified with the coordinate x~^ and a is a coordinate 
T-dual to z. 

There are 7 scalar fields rotated by an S0{7) symmetry. There is an additional 
scalar field x distinguished by the orbifold identification ^ and a U{N) gauge field with 
field strength F. Each of these adjoint- valued fields has canonical mass dimension 1 while 
a, T and L have length dimension 1. These are the natural dimension assignments from the 
perspective of gauge theory. Lastly, we note that as A 0, we recover standard matrix 
string theory [34-36]. 
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The singularity that develops in space-time as L — > is reflected classically in the matrix 
model by the appearance of a new flat directions at r = 0. These flat directions correspond 

to a fluctuations of the x' fields which are suppressed by the couplings, 



which vanish at L = r = 0. We will examine how this classical picture is modified quantum 
mechanically. 

This model can also be extended to describe string theory on the nuU-brane, or more 
generally M-theory on products of tori with the nuU-brane. This extension is straightfor- 
ward and involves promoting some of the x^ to gauge-fields while increasing the dimension 
of the field theory in the usual way [37] 



These extensions are already interesting for the case of string theory on the nuU-brane. 
Type IIA on the nuU-brane involves an extension of type IIB matrix string theory [35,38] 
in which time-dependent non-perturbative field theory effects should play an important 
role. On the other hand, type IIB on the null-brane involves a 3 -|- 1-dimensional time- 
dependent gauge theory generalizing [39,40]. The action of S-duality in this theory should 
be quite fascinating. We will not pursue these directions here. 

Lastly wc note that the M-thcory on the null-brane has supcrsclcction sectors specified 
by the quantized amount of graviton momentum on the z circle. This circle is large at 
infinity so these excitations are finite energy. Each sector with non-zero momentum can 
be studied in the matrix model by analyzing the quantum dynamics expanded around a 
classical background with non-zero electric fiux. 

2.2 Decoupling revisited and the Icirge N scaling 

The decoupling limit for the null-brane matrix model can be understood along the lines 
of the original argument given in [15]. Let us return to the orbifold perspective where 
we consider an array of DO-branes to make the connection between the arguments and 
determine how to take the large limit. 

The DO-brane action in string units has the schematic form 



{L' + lx'r'){D,xr, 



(7) 



X' 



(8) 
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The matrix theory hmit sends the dimensionless parameter gs to zero, holding fixed the 
dimensionless quantities 

i = gyh/is = Rt/el, X = g;^'^x/Es = x/ip. (10) 

The action expressed in terms of these parameters refers only to quantities in the resulting 
eleven- dimensional theory 



S= dt 



y — -j + [x, x] + fermions . (11) 



At these scales the interactions are of unit strength. 

We now come to the issue of scaling the parameters (L, A) of the null-brane orbifold. 
First of all, we should hold fixed the size of the neck in the geometry at r = in eleven- 
dimensional Planck units. The scaling of A is can be understood from the underlying DO- 
brane picture which is related to the 1+1 QFT appearing in (jHl) as a T-dual representation. 

In the DO-brane picture, one has groups of DO-branes in relative motion on different 
fundamental domains of the covering space of the orbifold. Strings stretching between 
image branes are represented by the Fourier modes of the 1+1 field theory. The DO-brane 
representation is, however, useful for understanding the matrix theory scaling limit of the 
null brane orbifold. The kinetic energy of relative DO-brane motion is excess energy above 
the BPS threshold 



Ey^ = E- N/R = ^{N/W+¥ - N/R = + ... (12) 

where the subleading terms in the last expression vanish in the scaling limit. The time scale 
(fnH) means the energy in the DO-brane dynamics should scale in order to keep Eym^sQs 
fixed. Equivalently, we are holding fixed 

ly = Eym ■ 2N{il/R) . (13) 

Now, DO-branes on neighboring images of the covering space have relative kinetic energy 

il{p-g-pf = {2n\ip/Rf. (14) 

In order to hold this quantity fixed, A must scale as A ~ gs''^. This is the scaling limit 
proposed in [20]. In particular, standard matrix theory arguments for the decoupling of 
both higher string modes and ten-dimensional gravitational physics may be directly carried 
over to the null-brane matrix model. 
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More precisely, decoupling holds as long as one takes the scaling limit first holding fixed 



(15) 



for fixed N. One can then contemplate further taking L/ip to zero, or to infinity, in 
various combinations. 

Consider now the large limit. In order to hold the invariant mass ()13p fixed as we 
scale N, the Yang- Mills energy Eym should scale as 1/A^. Since there are N DO-branes, 
the kinetic energy v'^/R of each must scale as l/N"^, or ~ 1/A^ {i.e. the velocity slows 
canonically as one boosts the frame by dialing A^). The relative velocity of DO-branes 
related by g'"' is 27mA; thus A ~ 1/A^. 

2.3 Fermion couplings from the quotient action 

To derive the fermion couplings, we can use the orbifold description of the nuU-brane. 
We will begin by studying the Euclidean instanton problem and then extend the result to 
DO-branes moving on the nuU-brane. The Euclidean instanton problem was also studied 



where are bosons with /i = 0, . . . , 9. The fermion, 0, is the dimensional reduction of a 
Majorana-Weyl spinor in ten dimensions. Both and G are in the adjoint representation 
of the gauge group. We will take the gauge group to be U{N). For the moment, we will 
take all fields to be dimensionless. 

We can express the action of the null-brane quotient group on the X^^ in the following 
way 



in [41]. 



The supersymmetric instanton action in fiat space is given by 



S = -Tt (- [X^, X,] [X^ - 07^ [X^ 6] 



(16) 



X exp(to')X + a 



(17) 



where 



UJ' 




(18) 



and 



LOx^ = —UJ+x = SttA, 



= 27rL, 



(19) 



with all other components vanishing. 
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The actions on the fermions can then be easily computed. The translations have no 
effect on a spinor representation, while the Lorentz transformation acts via 



e ^ exp (^luj^uY"^ e = (1 - 7rA7+") 6. 



(20) 



Note that half of the fermions, those that are annihilated by 7^, are invariant under the 
quotient action. This is equivalent to the statement that the background preserves one-half 
of the supersymmetries. 

2.3.1 Instantons on the quotient space 

We are now ready to study instantons on the quotient. This should really be viewed as a 
formal exercise on route to determining couplings for dynamical D-branes on the nuU-brane. 
So we view X^^ and ip as infinite matrices satisfying the constraints imposed by invariance 
under the quotient action 

m+k,n+k ^mn^ 



Xm+k,n+k — ^mn + 2vrA;AX^,, 
r — 

^m+k,n+k 



= X-„ + 27rA;AX„„ + 27r2A;^A2X+,, (21) 



^m+k,n+k ^mn ~l~ '^'^kLd^ymi 

Qm+k,n+k = (1 - 7rA;A7+'=) Qrnn- 

As usual, we Fourier transform this system (see [20] for this procedure applied to the bosons 
in this case). The matrix Qmn is replaced by the operator 



6{a - a') (22) 



where < a < 2tt. The function 6 {a) appearing in ()22j) is not Hermitian but this can be 
corrected by defining 

9 = 9+ '-^^^^9', (23) 

where the prime denotes a cx-derivative. Expressed in terms of 9, the operator G is Hermi- 
tian: 

e = ^ - ^7+^^' - y7''"^5<x- (24) 
Under a gauge transformation, 9 transforms according to the rule 

9 u9u^ + ^7+^ (u'9v) - M^M^') (25) 
10 



where u{a) is an element of U{N) for the case of N instantons. This is an unconventional 
gauge transformation rule. With L 7^ 0, we can define a natural gauge covariant fermion, 

^ = ^~+^7+'^,^~|. (26) 



4L 

If we also define the covariant derivative, Di = — iL^^z, then we see that 

e = ^-^7+-{Di,^} (27) 

which is manifestly Hermitian and covariant. 

We can combine this expression with the analogous operator expressions for the bosons 
found in [20] 

X- = X- + y {X, D,} - yDi ■ X+ ■ Di, 

Z = iLDi, (28) 
to compute the Yukawa coupling appearing in p6|) 

Tr (6)7^ [X^, 6]) = ^ j daTi (^LO^'DiO + 6 [-7-3;+ - 7+x" + 7^x + fx^ 

+T (^7+"^ [d. D^x+] + ^7+" [e, Dix'] - Di^7^ [e, x+] 

+^7^ {Di^, x+} + Di^7+ {0, x} - 9-1+ {D16, x}) 

\2 _ 

+— [Dih^ [Die,x+] -Die^+ [e,Dix+] 

o 

-e^+[D,e,D^x+])^. (29) 

We have dropped the hats on each 6 appearing in to avoid notational clutter. Some 
of the terms above are not manifestly Hermitian but this is not problematic since they can 
be made Hermitian by adding total derivatives. We will see momentarily that those terms 
vanish in the decoupling limit. 

In the abelian case where = 1, this Yukawa simplifies to the form 

e7„ [x^, 9] = — j da (iie-f'Die + — (^7+"^^Dix+ + e-f+"eDix' 

2n J \ 2 ' 

+2x+^7^Di0 - 2x^7+^1^)^ . (30) 
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It is worth noting that only couphngs involving two fermions appear in the action. No four 
fermion couplings are generated. This is natural since since we are considering a smooth 
quotient of flat space which has vanishing curvatures. 



2.3.2 The extension to quantum mechanics 

Let us recall the procedure which we followed in the boson sector to go from the instanton 
action to the Lorentzian theory describing DO-branes. We let our ten matrices be 
functions of r and we added a kinetic term D^X^D^X^ to the action. This term was 
reasonable because it was invariant under the ten-dimensional Poincare symmetry, if we 
assumed that r itself did not transform under this symmetry. Finally we gauge-fixed one 
of the coordinates 

This procedure has the advantage of being completely covariant until the final step, while 
also reproducing the correct couplings in all static cases. Note that r like a and L is 
currently dimensionless. We will restore all dimensionful factors below. 

Let us try to do the same with the fermions. We first note that in the quantum mechan- 
ics describing the DO-branes moving on the quotient space prior to Fourier transforming, 
the fermion kinetic term descends from the ambient fiat space expression. All of the compli- 
cations arising from the Fourier transform are captured by the Yukawa coupling (j^^ . We 
can then apply the decoupling limit to (j^^ which amounts to scaling A ~ e and dropping 
terms that vanish as e ^ 0. Let us call the 1 -|- 1-dimensional fermion ip. The proposed 
1 -|- 1-dimensional action after decoupling is given by 

S = J j drdaTi (^{Dox'f + {D^xf - L\Dix'f - L\Dixf - V2\{x - tDqx)F 

1 1 9 ?7~A — — 

+ + -AV2)F2 + - [x\ x'] +{[x, x'] + ^D^x'f + z^70DoV' + tL^^'D^^: 

2 2 a/2 

+4j [7"x + j'x\ 4j] - '-^D.^Yt/j) (32) 



V2 

where i = 1, . . . , 7 and we have set = 1. In the abelian case, the action simplifies to 
become 

S = - I drda ({dox'f + {doxf - {L^ + ]-X^r^){d,x'f - L^{dixf - V2X{x - Td^x)F 
9s J \ 2 

+ {L' + ^X't')F^ + i^p-f'^doiJ + i^LY + ^7ldii?j . (33) 
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Note that the matrix + (rA/-\/2)7^ squares to give (L^ + ^A^r^) times the identity 
matrix. It should also be possible to derive the fermion couplings from the DBI action 
taking carefully account of the S-field and dilaton along the lines described in [42-44]. 

3 The Leading Quantum Corrections 

We now turn to quantum corrections to the effective dynamics. The moduli space for the 
1 + 1-dimensional theory defined by (jH^ is parametrized by a choice of vacuum expectation 
value for (x*, x) and a choice of Wilson line in the a direction. We will put aside the choice 
of Wilson line for the moment. 

For example, for gauge group SU{2) we parametrize the moduli space by the Spin{7) 
and gauge-invariant combinations 

(Tr(xO^) = ^, (Tr(x2)) = ^. (34) 

The choice of a U{2) gauge group describes the interactions of two gravitons moving on the 
nuU-brane space-time. Restricting to SU{2) factors out the center of mass motion. 

Now a glance at (jH^ shows that the gs plays the role of h. We will compute the static 
potential energy as a function of (6, c) to leading order in gs- To do so, we employ the 
background field method. 

3.1 Gauge- fixing the action 

The first step needed to compute the quantum corrections to the action is to fix the gauge 
symmetry, and determine the ghost content of the theory. We need to implement a back- 
ground field gauge condition. The usual gauge-fixing condition used in [45] for DO-branes 
in flat space is 

G = -iAo + v,u[X^,X'']=0 (35) 

where {Aq,X'^) are fluctuating fields while are background fields. This choice of gauge 
fixing term ensures that the resulting quantum effective action is invariant under gauge 
transformations of the background fields. 

Into this expression, we insert the operators appearing in (j2Hl) after taking the decoupling 
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limit. The resulting gauge-fixing term takes the form 

G = -iAo + iAA[ + -^x' + [xl x'] + [xb, x] + ^ [xb, Ai] 

V 2 V 2 

+^ [{A,)b, x] + [iAo)b, Ao] + A Pi),, A,] , (36) 

where A = + |A^r^. Here a dot represents differentiation with respect to r, while a 
prime indicates differentiation with respect to a. 

The ghost action is obtained from the variation of G with respect to a gauge transfor- 
mation acting on the fluctuating fields. Explicitly for the case at hand, the infinitessimal 
transformations with gauge parameter u are 

6x' = i [0, x' + , 

5x = i[0,x + Xb] , (37) 

Ml = 0' + z[0,Ai + (Ai)b], 

6A0 = + i[0,Ao + {Ao)b]. 
Finally one also adds a gauge-fixing term to the action (p?^ given by 

<S'gaugc-fixing = — drdcr Tr (G^) . (38) 
9s J 

3.2 Warm-up: the abelian case 

Our real intention is to compute the effective potential for the moduli fields b and c. How- 
ever, as a warm-up exercise let us consider the abelian case where we expand around the 
constant background configuration 

Xb = xo, xi = xl {Ao)b = 0, {Ai)b = 0. (39) 

In this case we use the gauge-fixing term, 

G = -iAo + iAA[ + —x', (40) 

v2 

which varies under a gauge transformation in the following way: 

6G = i{-0 + A0") . (41) 

This leads to the ghost action 

5'ghost = - / drda C* (-52 + Adf) C. (42) 
9 s J 
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This, in part, explains the time-dependence appearing in the gauge-fixing term (jHUj) . The 
ghosts, hke the gauge fields, see a time-dependent a circle. 

To obtain the other pieces of the action, we replace by + Xq in and add the 
gauge-fixing term. The result is 

S = 1- j drda I {x'Y + x^-Al + 2AAoA[ + V2t\Aox' - A {x"Y - A (x'f - A^ {A[f 
-V2r\AA[x' - V2\ {x + Xq - rx) F + AF^ + i^j-i^i) + i^) (^L^ + ^^""^ (^3) 

where 

F = Ai- A',. (44) 

In this abelian case, the coupling dependence appearing in is irrelevant so we will drop 
the 1/gs factor. 

The bosonic part of this action has non-standard kinetic terms. To correct this, we must 
perform a field redefinition. Define 

x = x + ^Ai. (45) 
v2 

This definition does not respect gauge-invariance but since we have already gauge-fixed, 
this is not a problem. With these redefinitions, the bosonic action simplifies to 

S = J drda I {x'Y + S:^ + Al-Al-A [(a;'')' + {x'f + {A[f - {A'^fj + 2V2XxA'^'j . (46) 

We decompose the fields into their Fourier modes labeled by n. At each level, we find 
eight bosons of "mass" = —n'^{L'^ + A^r^/2) with conventional kinetic terms. The 
operator governing quadratic fiuctuations is 

H=-d'^-n'(L' + ^] . (47) 



2 

Integrating over the quadratic fiuctuations for each such boson generates a factor, 

exp (^i J drda Ks^ = det-'^^{H) = det'^/^ |^_^2 _ ^2(^2 ^ . (48) 

To express this as an effective potential, we use a proper time representation for the deter- 
minant, 

det-'/\H) = exp / ^ / ^e--^^-\r, r)) , (49) 
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where t is the Lorentzian proper time, and an ie is inserted for convergence. The integral 
over 0" is not part of the trace but can be inserted harmlessly in this case because a does 
not appear in the kernel. 

The only information we require is the propagator for H, given by a continuation of 
Mehler's formula for the simple harmonic oscillator 

e-^*^(n,r2) = e-^^^*f ^ ^ V^' (50) 



X exp 



'Ut^ + r|) cosh( V2nAt) - nrs 



.y2sinh(v/2nAt) ^ ^ 

The large t behavior of ()49p corresponds to the infrared contribution to the potential. 
From ()50|) . we see that this contribution is highly suppressed by both the n^L^ mass term, 
the time-dependent nr^ term, and the sinh(\/2nAt) prefactor. This is what we expect for 
a massive particle with a growing time-dependent mass. 

Now there are two additional bosons x and Aq. This is a coupled system that is more 
subtle to analyze. Path-integrating over x and Aq results in the determinant, 

det-'/\H^ + q^) (51) 

where = 2r;,^A^. Note that the spectrum of + is well-behaved because H is Hermi- 
tian, and gapped because of q. There is a convenient way to represent the determinant 
using the product, 

det-^^^{H + iq) X det-^/'^{H - iq). (52) 

The reason this is convenient is that we can directly use the kernel ()50|) to represent the 
effective potential. This results in complex masses, 

= n'^{L^ + ^-) ± V2in\, (53) 

which really require different proper time representations ()49|) (differing in the sign of t) 
depending on the sign of the imaginary part of the mass. We will avoid this issue and 
directly compute the determinant (jSH). 

The complex ghosts with action simply cancel the contribution of 2 of the 8 bosons 
with mass = n^(L^ + r^A^/2). We are left with the fermions. To find the fermion mass 
spectrum, it is convenient to introduce complex masses much as we could have done in the 
(x, Aq) system. The path-integral over the real fermions gives the Pfaffian of the operator, 

-tdr + ti'' (^LY + ^r'^d^. (54) 
16 



This is a Hermitian operator with real eigenvalues. It is more convenient to rewrite this 



Pfaffian as y det(^) in terms of the Dirac operator. The square of the Dirac operator unlike 
the square of (jH^ takes a nice form: 

^7°5. + ^ {li^ + ^7") d}j = dl - Adl - ^7°7"5.- (55) 

The eigenvalues of the matrix n^A — [inX/ v^)7°7^ are complex: n^A + in\/\/2 and n^A — 
in\/\/2. The 16 real fermions therefore contribute for each n a factor 



det^ 

which we express as 



dl + n^A - in\/V2 



X det' 



det^ [H^ + gV4] . 
Collecting together these results, we need to compute 



det-i/2 



exp Tr 



2\og{H' + q'/A)--\og{H\H^ + q^)) 



(56) 
(57) 

(58) 



L(i/2 + g2/4)"J 

to determine the potential. We will expand the logarithms formally in inverse powers of H 



3 b q^ 63 q^ 51 q 



10 



+ 



+ .... (59) 



2-ov- V- . . ^g^4 32 i76 512/78 5i2ifio 
Now we can represent inverse powers of H using a Lorentzian proper time formalism 



This amount to the replacement 



(n- 1)! 



1 {itY 



t {n - 1)! 

in the series ()59|1 . On substituting, we find the series 
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3^(M* + ^M« + ^(rf 



;(?«) 



10 



sinh^ ( — 
t \ 4 



(60) 



(61) 



(62) 



1935360 

This results in the following Lorentzian proper time expression for the effective potential 

JrfiL'^t 

o biiili 1 — — I r 

V 9./9 / 

n=0 



eS 



^ /^VSsinh^ (^) 
27r y t \2V2J 



X 



inX 



2^^^/2smh{^/2nXt) 
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1/2 



exp 



inXr 



tanh(^)]. (63) 



Note that this potential ()63|) has an apparent large t infrared divergence which we will 
revisit after rotating to Euclidean space. 

There are a few points worth noting. Let us rewrite the 8 sinh'' j term in the form 

^ ^8 - 2 + 2 cosh{V2nXt) - 8 cosh (^^) . (64) 



This is the collection of prefactors we would have obtained had we directly represented 
and ()56j) using the same proper time representation regardless of the sign of the imaginary 
part of the mass. If we restrict to n = 0, we see that there are effectively 8 + 2 — 2 net 
bosonic determinants which are precisely canceled by the fermionic determinants. This 
reflects the underlying supersymmetry present in the n = sector. 

Similarly, the leading UV divergence also cancels in the effective potential. This comes 
about because the cosh prefactors that distinguish the fermion and {x, Aq) effective potential 
terms do not change the leading t — * behavior of the determinants; the underlying 
supersymmetry still kills the divergence characteristic of 1 + 1-dimensional field theory 
in the proper time integral in ()49j) . 

We can now rotate to both Euclidean proper time and Euclidean time sending 

t —>■ it, X —>■ —iX, T — »• ir. (65) 

The reason we need to rotate A is to ensure that the sinh factors in continue in a 
sensible way to Euclidean space. One finds 

1/2 



n=0 



X exp 



2^2/ y27rv/2sinh (v/2s) 

nAr^ . f s 



tanh 



(66) 



V2 \V2, 
where we have defined s = nXt. 

We see that either nr'^X ^ 1 or nlP'/X ^ 1 forces s -C 1; in other words, nA/X ^ 1 
implies only s ^ 1 contributes to the integral. In this regime, we have 

ds . „A„/\ /'nX\ 1/2 



2vrK. ~ /^,4y^-„A./A/!^ 
J s ^ Vtts 



- ^C(3)r(7/2)-^ (67) 

This result is easily understood from the DO-brane picture as arising from the one-loop 
y4^^7 interaction of matrix theory, summed over images under the orbifold group. 
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When s becomes sufficiently large, the effective potential (jUBj) can have an infrared 
divergence. This occurs roughly when 

nL^ < X/V2. (68) 

From PUI) . we see that the right hand side is the coupling constant in the abelian theory 
mixing the modes of {x,Aq), so the divergence occurs when the energy of the lightest 
Kaluza-Klein mode is of order the scale set by the coupling constant. Strong coupling 
physics appears in the low energy regime of the matrix model, where the dimensionful 
Yang-Mills coupling becomes effectively large. The infrared divergence appears to herald 
the onset of strong coupling, and the breakdown of perturbation theory. 

Given that the ultraviolet contribution to the potential collapses to the f ^/ r'^ interaction, 
it seems quite possible that the non-renormalization theorems for this and higher velocity 
interactions [46] can be extended to this potential. 

The result (p7j) yields a time-dependent contribution to the energy of the center-of-mass 
dynamics of the system that depends on the orbifold parameters (L, A) but is independent 
of the center of mass position and velocity. This energy indicates the time-dependent scale 
of supersymmetry breaking in the system. 

3.3 Particle production 

A time-dependent background generically results in particle production. While there is 
no particle production in the null-brane vacuum because of the existence of a null killing 
vector, there is particle production in the matrix description. This is because the matrix 
theory describes objects carrying longitudinal momentum which preserves a different set of 
supersymmetries than the null-brane background. 

The linearized dynamics of the bosonic fields is solved by parabolic cylinder functions 
(our discussion here follows [47]), 

^^ = «(^^---H^"^^) + ^(M^^Vil-^^-)' 

where the Bogolubov coefficients {a,P) satisfy |ap — = 1, and 
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(69) 



(70) 



We use the values u = n X/V2 and 



'L"^ relevant to the dynamics of the ra*^ mode of 



X*. Note that fi is simply the dimensionless ratio of the energy cost and the coupling of a 
given mode, see equation 

The parabolic cylinder functions have the asymptotics 



D 



2/1- 



.1 I e 4 



2UJT 



/2uiT 



T 



+00 



g i(fr2+;.Iog(v/^|r|)) 



(71) 



where 7 is the phase 



7 



r(|-^/i) 



r 



-00 



4 log 



+ ;Ulog2. 



(72) 



The Bogolubov coefficient f3 expressing the probability amplitude for particle production 
may be read from (|7T|): 



-TTfJ, 



exp 



TinL^ 
V2X 



(73) 



Thus the mode production rate becomes of order one as the effective coupling of a mode 
becomes of order one. This feature is related to the onset of infrared divergences in the 
one-loop determinants calculated in the previous subsection; the infrared singularities and 
the copious mode production are both signals of the breakdown of the weak-coupling per- 
turbative expansion. 

The above calculation of mode production is nicely reproduced from the underlying 
picture of DO-branes on the nuU-brane orbifold. Consider a DO-brane and its image under 
the orbifold group They undergo a scattering process which at leading order in string 
perturbation theory was calculated in [48]. The real part of this annulus amplitude can be 
interpreted as the eikonal phase shift of scattering D-branes (see for instance [48,49]); the 
imaginary part of the phase shift gives the probability of vacuum decay. Alternatively, the 
imaginary part is the pair production probability of open strings stretching between the 
branes via the optical theorem. These stretched strings are the Fourier modes in the T-dual 
field theory, and so the production amphtude of these stretched strings should agree with 
the field theoretic calculation above. 

In the matrix theory limit, the annulus amplitude becomes the one- loop amplitude in 
the field theory ()32|1 . leading to the collection of determinants evaluated in the previous 
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subsection (integrated over r). The result of [48,49] for the scattering phase shift of a 
DO-brane and its image under (7" is 



5n = ~- 



tanh" 



1 rd^-snLVX. 



sm s 



f4cos(s/y2) - cos(y2s) - 3 
i/y2) V ' ' ' 



v^mL^ + A 



-■ 2 



(74) 



(75) 



The result fl74|) is simply the local expression for the determinants in r, given in equa- 
tion ()66p. integrated over r and continued back in A. The full phase shift sums this result 
for 6n over the images labelled by n, which is the same as the sum over modes in the T-dual 
field theory. 

How is this decay probability related to the Bogolubov coefficient calculated above? 
The decay probability is the overlap of the in- and out-vacua, | out(0|0)inp. The Bogolubov 
transformation of (uncharged) bosonic creation/annihilation operators 



flout = actin + 

(with |ap — = 1) allows one to write the out vacuum as 



|0) 



out 



a 



exp 



. . 



|0)i 



and thus 



-log(|out(0|0)i,|') = log|a| = |log[l + 
The analogous calculation for fermions yields, 

-llog[l-\/3\% 

because in this case the Bogolubov coefficients are related by |ap -l- 



(76) 



(77) 



(7J 



(79) 



1. Finally, 



log 



1 - 



logtanh(log|/?|"^) 



^0) 



connects the decay probability to (fTSj) . Note also that we can understand the real part of 
the phase shift 6 as the phase 7, appearing in (|7^ . accumulated in propagating over the 
inverted oscillator barrier (taking into account the frequency shifts involved for the various 
bosons and fermions). 
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This relation between the two results (fTSjl and (f7S|) yields /5 = e~'"'^ = e"'^''^/^''. The 
group element g"' of the orbifold generates an image DO-brane whose impact parameter is 
b = nL and relative velocity v = n\/ 

We see in the result fl73|l the effect of the order of limits, and the various scales. If we 
send L — > 0^ at fixed A, mode production is of order one for all the A^^ modes of the system. 
The dynamics after passing through the neck of the geometry involves energies that grow 
with iV and will not have a uniform large limit. On the other hand, this is a DLCQ 
artifact. Recall from our discussion in section IT^ that A ~ because it is a velocity in 
the boosted frame of the matrix model. 

Taking the large N limit first, particle production is an effect suppressed by e~'^^ for 
some constant c, and thus totally unimportant. Note, however, that strings that thread 
through k N DO-branes have an action cost that is not suppressed parametrically in 
A^; these strings will be produced in the passage through the neck, however this effect is a 
collective multi-particle excitation from the point of view of the perturbative matrix model, 
not governed directly by the production rate ()73|) . It is an open question whether such 
collective excitations are well-behaved in the large N limit. The answer to this question is 
central to the issue of whether matrix theory gives a well-defined answer to what happens 
at the cosmological singularity at small L and large N . 

3.4 The non-abelian theory 

Now let us to the case of a broken U{2) gauge group. We expand around the vacuum 
configuration ()34|) generalized to permit velocities with the explicit choice of background 
fields 



Ad ^2 



where a" are the Pauli matrices. The unusual form of x\, is required to satisfy the classical 
equations of motion which include mixing of x and F . However, if we again make a change 
of variables to x = x -|- \tA\I -\/2, then we have the simpler setup 

x, = ^^a\ 7 = 7 + ^- (82) 

Since the dependence on the relative velocities 7, and 5 follows from Galilean invariance, 
we suppress the velocity dependence for the moment but will restore it later. We then 
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consider the fluctuations around this background, 

1 



+ etc. 



^3) 



As usual, we will we ignore the U{1) center of mass physics and focus on the interacting 
non-abelian theory. 

We will repeat the procedure used in the abelian case except we will work directly with 
complex masses. We flnd a collection of particles with masses conveniently parametrized 
in terms of 



± 



L'^{d±n) 



b'+ c± 



nXr 

7f 



^4) 



For each n, there are 8 massive bosons with masses, = and 8 with masses = r^. 
There is also one complex ghost with = and one with = r^. The shift in r 
induced by c in ()84j) is physically signiflcant when we allow the moduli flelds to vary slowly 
as functions of space-time. 

In addition to these particles, there are again a pair of real bosons with = r 



± \/2inX. The complex masses appear for exactly the 



and a pair of bosons with = 
same reasons as in the abelian case. 

Finally, the fermions also split into two groups with the first 16 real fermions generating 
a factor of 



det^ 



92 + 



\/V2 



X det' 



d^ + rl + in\/V2 



while the second 16 real fermions generate the determinants 



det^ 



92 + r? 



tn 



\/V2 



X det' 



dl + r'^ + in. 



The potential therefore splits into the sum of two contributions. 



i5) 



i7) 



where 



eflf 



^ J y f^(^8-2 + 2cosh(y2nAt)-8cosh(^^^^i^±(t,^±Ar,n), 



dt 



n=0 

00 



n=0 



nXt 
2V2 



n 



using Euclidean proper time, and 



K±(t, T, n) = e 



-{{n±dfL'^+b'^)t 



27rv/2sinh(y2nAt) 
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1/2 



exp 



■ tanh 



nXt 

7f 



^9) 



Note that from the discussion in the abehan case, we see that the potential is UV finite 
{i.e. as t — >■ 0). There is an infrared divergence when 

{n ± dfL^ + b^< n\/V2 . (90) 

As in the abehan case, this appears to signal the breakdown of perturbation theory. 

Once again, for nr'^/v ^ 1 the dominant contribution is from the UV region t 0. 
As in the abelian case, we can understand the effective potential as due to the sum over 
orbifold group images of v^/r'^, where 

rl = {b + PTf + {c + ^T±^T)^ + {d + 6T±nfL\ (91) 

and we have restored the full dependence on the velocities. We have been a little cavalier 
about the sign of the potential. However, from the DO-brane picture, we can see that a 
potential that originates from the v^/r'^ interaction will have the same sign as the kinetic 
terms in the action [50]; the potential is therefore attractive. 

The potential is decaying rapidly with |r|. This implies that the fiat directions are 
restored as |r| oo and we recover the conventional picture of gravity in matrix theory, at 
least in terms of the structure of the potential. We can get some feel for the potential by 
setting c = d = along with the velocities = 7 = 5 = and A = 1. The corresponding 
static potential as a function of r appears in figure El The static potential as a function of 
impact parameter for various choices of L appears in figure El 

3.5 Final comments 

There are a few obvious generalizations of the above analysis. As we described briefiy in 
section \'2.1\ one can compactify additional dimensions. For instance, a toroidal compacti- 
fication in p additional directions will extend the 1+1-dimensional field theory to a theory 
on T^+i, one of whose cycles has time-dependent proper size A. The most obviously inter- 
esting cases are p = 1 and p = 2 where we expect novel field theory phenomena to occur: 
for p = 1, time-dependent instanton effects while for p = 2 we expect to find a version of 
S- duality. 

However, there are other ways to compactify. We could have all the cycles undergo a 
bounce like the bounce in the z-direction by combining the translational identification with 
a boost as in ((H); then each cycle has a time-dependent size {Lf + |Afr^)^/^. 
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An object passing through the neck of such a bouncing generates a singularity as we 
send the Li to zero. Again, for small but nonzero Lj, what initially forms is a black p-brane 
filling the torus. As the torus expands after the bounce, this p-brane will experience a 
Gregory-Laflamme instability and break up into a collection of black holes. This result is 
reminiscent of the proposal of [51] for a cosmology whose initial state is a dense gas of black 
holes. Of course, in the present case the gas of black holes is anisotropic, filling only those 
directions that are compactified; we do not currently have a non-perturbative description 
of the situation when all spatial directions are compact. 
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